a^ 



Ph. 



en 
> 

a^ 
in 
^, 

i> 

o 



X 



Matrix Superpotential Linear in Variable 

Parameters 



Yuri Karadzhov 

Institute of Mathematics, National Academy of Sciences of Ukraine, 
3 Tereshchenkivs'ka Street, Kyiv-4, Ukraine, 01601 

o 

(N 

C^ ; Abstract 

^^ , The paper presents the classification of inatrix valued superpotentials corre- 

sponding to shape invariant systems of Schrodinger equations. All inequivalent 
irreducible matrix superpotentials realized by matrices of arbitrary dimension 
with linear dependence on variable parameter are presented explicitly. 



^E-mail: yuri . karadzhovOgmail . com 



1 Introduction 

Supersymmetric quantum mechanic presents a powerful and elegant tool for obtaining 
explicit solutions of quantum mechanics problems described by Schrodinger equations 
m . Invented by Gendenstein [2] , property of discrete reparametrization of potentials, 
known as shape-invar iance, helps to determine whether eigenvalues of Hamiltonians 
can be calculated by algebraic methods. 

Though shape-invariant potentials do not exhaust the full class of potentials of 
solvable Schrodinger equations \13\, it was interesting to find new integrable models. 
Previously several attempts of describing the class of shape-invariant potentials were 
made. In paper of Cooper et al [3] a wide class of scalar shape-invariant potentials 
was described. 

In the world of matrix superpotentials there is still room for further research. An 
attractive example of a matrix problem which admits a shape invariant supersymmet- 
ric formulation was discovered by Pron'ko and Stroganov [4J, who studied a motion 
of a neutral non-relativistic fermion which interacts anomalously with the magnetic 
field generated by a thin current carrying wire. The supersymmetric approach to 
the Pron'ko-Stroganov problem was applied in papers [5], [6], [7]. Particular cases 
of matrix potentials were discussed in papers [S] , [S] , [ID] , [S] • Matrix superpoten- 
tials appear also in some supersymmetric systems related to the crystalline structures 
in Gross-Neveu model [IT], [18]. Two-dimensional matrix superpotentials, including 
shape invariant, were studied in papers [15], [16]. In paper pj^ Fukui considers a 
certain class of shape invariant potentials, which however was ad hoc restricted to 
2x2 matrices with superpotentials linearly dependent on the variable parameter. 

A systematic study of the problem was carried out in our recent paper [12] where 
we presented a complete description of irreducible matrix potentials, which include 
terms linear in variable parameter, and are proportional to the unit matrix. 

In the present paper the classification is continued and Fukui's results are gener- 
alized. All shape-invariant matrix potentials of an arbitrary dimension with super- 
potentials linearly dependent on the variable parameter were found. 

2 Shape- invariant potentials 

Let's start with a spectral problem 

Hki/J = EkiJ, (1) 

where Hk is a Hamiltonian with a matrix potential, Ej^ and -0 are its eigenval- 
ues and eigenfunctions correspondingly, moreover ip is n-component spinor. In the 
Schrodinger equation Hamiltonian has the form 

Hk = -^ + V,{x), (2) 



where Vk{x) is a ra-dimensional matrix potential depending on variable x and param- 
eter k. 

Suppose that Hamiltonian accepts factorization 

Hk = alak, (3) 

then its superpartner is defined as follows 

H+ = a,al (4) 

The most common representation of operators ak and a|. has the form 

d d 

ak = Ak{x)— + Bk{x), al = -—Al{x)+Bl{x), (5) 

where Ak{x), Bk{x) are matrices depending on x and ^4^(2;), BI(x) are hermitian con- 
jugate to them. 

Substituting this representation into equation ([3]) we obtain the equation 

H, = -A\A,— + {BlA, - AiB, - (AIA,)')— + bIb, - (4i?,)', (6) 

which is supposed to be Schrodinger equation of the form (I2j). It leads to the following 
conditions 



(7) 



In terms of new variable Wk{x) = Al.(x)Bk{x) this condition take the form 

The same result can be obtained analogously with simpler representation of op 
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AlAk = I, 




BlA, - AlB, - {AlA,y : 


= 


BlBk - (AlBk)' = Vk. 





9 , (8) 

V, = Wl -Wi ^ ' 



erators ak and a. 



ak = ^ + Wk{x), ai = -— + Wk{x). (9) 

Wk{x) is called a matrix superpotential. 

As Wk{x) is hermitian, then the corresponding potential and its superpartner 
y+(x), i.e. 

dWf, o , dWk n , , 



are hermitian too. 

The goal is to find such superpotentials which generate shape-invariant Hamilto- 
nians 

Ht = HF^ + C,, (11) 

where F^, Ck are scalar functions of k. In terms of superpotential, the last condition 
has the form 

Wl + Wi = Wl - W'p^ + C (12) 

It is sufficient to search for irreducible superpotentials, which means matrix W^ 
can't be transformed to block-diagonal form with a constant unitary transformation. 
Note that ii Fk = k then equation ([12]) transforms into 



W'k = \Ck (13) 

and it follows that 

Wk = \ckx + nk, (14) 

where Vtk is hermitian matrix depending on k which can be diagonalized. So in this 
case superpotential is completely reducible to the one-dimensional shifted oscillators. 
In case Fk ^ k it makes sense to restrict ourselves to unit shifts 

Fk = k + l, (15) 

see discussion section for details. 

In the following sections the classification of superpotentials linearly dependent 
on parameter k is presented. 



3 The determining equations 

Let us consider the superpotential of the form 

Wk = kQ + P, (16) 

where P and Q are n x n hermitian matrices dependent on x. 

Suppose that Q is not proportional to the unit matrix because in this case the 
superpotential is reducible. 

Substituting fITB]) into flT^ and taking into account equation flT^ we obtain the 
equation 

{2k + l)(g2 - Q') + {Q, P} - 2P' = Ck, (17) 



where Q' = ^^P' = ^ and {Q, P} = QP + PQ is an anticommutator of matrices 
Q and P. After variable separation the last equation transforms into a system of 
determining equations 

Q' = Q^ + V, (18) 

P' = \{Q.P}-i^. (19) 

Cfe = (2A; + l)z/ + 2/x, (20) 

where /i, ;/ G M are arbitrary constants. For scalar values unit matrix / is omitted 
throughout the paper and written, for instance, as /i instead of /i/. 

4 Solving the determining equations 

In this section a solution for determining system of arbitrary dimension n is presented. 
Let us start with equation flTSl) . 
Introducing new variable 

M = Q-Lp, (21) 

where scalar function ip is defined in the following way 

Atan(Ax + 7), z/ = A^ > 

-Atanh(Ax + 7), i/ = -A^ < (22) 



V' 






and 7 G M is chosen so that M is not a singular matrix, we can transform equation 
([IBD into 

M-^M'M-^ = 1 + 2ipM-\ (23) 

Note that 

(M-^y = -M-^M'M-\ 
So if we set A^ = M~^ then we obtain the following linear equation: 

A^' = -/ - 2ipN. (24) 

Its general solution has the form 

N = -p{x)I + e{x)C, (25) 



where p{x) and 6{x) are scalar real valued functions and C is an arbitrary constant 
matrix. 



P 



^sin(2(Ax + 7)), z/ = A^ > 
^sinh(2(Aa; + 7)), z/ = -A^ < 
X + 7, 1^ = 

cos^(Ax + 7), z/ = A^ > 
cosh^(Ax + 7), z/ = -A^ < 



(26) 



(x + 7)^ 



z/ 







Q is liermitian, so A^ is liermitian too. Moreover, if unitary matrix U is such 
that UNW is a diagonal matrix, then UQW is also diagonal. As p and 6 are scalar 
functions, then diagonalization transformation U should not depend on x. So Q can 
be diagonalized 



Q = diag {§!,..., g„} 
and equation f lTH]) transforms into system 



Qi = (i + T^, z = l...n, 



which has following solutions 

gj = Atan(Ax + 7j), i = l...n^ z/ = A^ > 0; 
— Atanh(Aa; + 7j), z = 
— Acoth(Ax + 7i), i = 
. ±A, i = 

Z = 1...7TT, 

m + l...n 



(27) 



(28) 



Qi 



Qi 



1 



0, i 



l...m 

l + l...n 
, z/ = 0, 



-A^ < 0; 



(29) 



where 7^ G M, i = l...n are integration constants. In cases i^ < and u = matrix 
Q consists of blocks of size m,l — m + l,n — l + l and m, n — m + 1 correspondingly. 
Some of them may or may not appear, i.e. have zero size. 

So we define matrix Q up to constant unitary transformation U. Now let us 
consider linear equation (TT9|) for P = {pij} which can be solved element-wise: 



Pa = -tan(Ax + 7i) + </?« sec(Ax + 7^), i = l...n 

Pij = (pij\/sec{Xx + 7j) sec(Ax + 7^), i = l...n,j = 1... 



(30) 



If z/ = -A2 < 



Pa 



-^ tanh(Ax + 7i) + Lpu sech(Ax + 7^ 



l...m 



-^ coth(Ax + 7j) + (^ii csch(Ax + 7i), i = m + 1...1 



_ ±^ + V9jjexp(±Ax), z = / + l.-.n 

(^jjA/sech(Ax + 7j) sech(Ax + 7j), « = l...m,j 



l...m 
m + 1...1 
l + l...n 



(31) 



V'jj V'sech(Aa; + %) csch(A2; + 7j), « = l...m,j 

ipij ^ysech.{Xx + 7j) exp(±Ax), z = l...m,j 

Pij = ipij ^ycsch{Xx + 7j) csch(Ax + 7j), i = m + 1..J, j = m + 1.../ 

y?ijA/csch(Ax + 7j) exp(±Ax), z = m + 1..J, j = / + l...n 

ipijexp{±\x), (*) z = / + l...n, J = / + l...n 

(/jy, (*) i = I + l...nj = I + l...n 

Use formula (*) if corresponding diagonal entries qi and Qj have the same sign and 
formula (*) otherwise. 

• If z/ = 



Pa 



Pij 



-ix+2'ri) 



-fix + (fa, 

^(x+7i)(a:+7j)' 






z = l...m 

z = m + l...n 
i = l...m,j = l...m 
i = l...m,j = m + l...n 
i = m + l...n,j = 171 + l...n 



(32) 



where (fji = Tp~J G C are integration constants. Numbers m and / in the above intervals 
correspond to the numbers defined in (12^ . Matrix P can be further simplified with 
unitary transformation, that reduces the number of non-zero entries in similar blocks. 
Obtained superpotentials are irreducible if there are enough non-zero entries tpji. 

5 Solving spectral problem 

Consider the Schrodinger equation 

Hki/J = {Hk + Ck)i^ = Ekij, (33) 

where Ck vanishes with constant multiplied by unit matrix in the Hamiltonian Hk- It 
follows from shape-invariant condition (TT2!) that 



Ck — Cfc+i — Ck- 



(34) 



Since all considered Hamiltonians are shape-invariant, equation fl5^ can be solved 
using the standard SSQM technique. An algorithm for constructing exact solutions 
of supersymmetric shape-invariant Schrodinger equations can be found in [3] . 

• Ground state '0°(x) is proportional to the square-integrable solution of the first 
order equation 

ak^l^lix) = (^ + W,\ ^°(x) = 0. (35) 

Function ip^ solves equation fl33l) with eigenvalues 

El = -Ck. (36) 

• Solution which corresponds to the n*'^ excited state ipki^) ^^^ ^e represented 
as 

C(a;) = 44+1 ■ ■ ■ (^i+n^ii^k+nix) (37) 

and must be a square-integrable function of x. The corresponding eigenvalue is 

-E'fc = —Ck+n- (38) 

It follows from equations ([20]), (El, dSSD, dMD that 

E'^ = El~ 2nfi - (n^ + 2kn)u. (39) 

At the end of the section, an example of exactly integrable problem is presented. Let 
Q and P, defined by formulas fl29]) . f lST]) . be 



i n\ I ^ - — -^ 

Q={ n^ n)^ P=\' 2^ ,.f 1> (40) 



OJ' I -^ lix 



where /i, (p are real constants, ^ > 0. Then corresponding Schrodinger equation looks 
hke 

and Cfc is equal to yU. 

The square integrable ground state solution "0° = I to ) has the form 

e = Cl6 + C2^2 ^ ' 



for non integer A; > 0, where 

01 = x^+le-'^HBik, 0, -2 - A;, ^, ^x), 

6 = f x^+ie-^//l,(A:, 0, -2 - fc, ^, f x)+ (43) 



2^^ ■- -'-'_BV"'?"i '' "■! ^1 2 



,2 



and 



^a;'=+V^ifB(fc,0,-2-fc,^,fx) 



62 = x5e-^//B(-A:, 0, -2 - A;, ^, ^x) 



6 = f xe-^//^(-^, 0, -2 - ^, ^, f a;)+ (44) 

-^e-^if,(-fc,0,-2-/c,^,fx). 

Here iif^ denotes the Heun Biconfluent function. For integer k solutions (143 p and 
fj4i|) are identical and second linearly independent solution of equation fl35l) should 
be used instead of fH4|) . It is however cumbersome and is omitted in the paper. 
Note that constants Ci and C2 should be chosen so, that 

00 

(V'°)V°^3; = 1. (45) 

■00 

The spectrum of the problem is described by formula 

E" = -(2n + l)/i. (46) 

6 Discussion 

In this section the possibility of transforming unknown function Fk in the shape- 
invariance condition flTTl) to unit shift is discussed. 

During the computation it was assumed that F^ = A; + 1, let us show it is a 
reasonable restriction. Under invertible transformation of variables 

k -^ a{k) (47) 

the function Fk = F{k) changes by a similar transformation 

F{k) -^ aFa-\k) = a{F{a-\k))). (48) 

Searching for such transformation, that would change function F[k) to unit shift, we 
get the equation 

F{a'\k)) = a-\k + l). (49) 

The above equation is known as Abel functional equation. The results concerning 
the solution of this equation were obtained in papers [19], [20], [21], [22] . 

Let X be R or R+. It is proved that 



(CI) if F : X ^ W is an injective function such that for every compact set K G X 
there exists p eN such that Vn, m G N°, \n — m\ > p : 

F'^iK) n F'"(K) = 

then there exists a solution to the Abel functional equation. 

So if the above condition is fulfilled, F^ can be transformed into unit shift. 

7 Conclusion 

Restricting the class of superpotentials to a simple form ( 1T^ . it was possible to find 
new exactly solvable systems of Schrodinger equations. As corresponding Hamiltoni- 
ans satisfy shape-invariant condition flT2|) . the systems can be solved using a standard 
SSQM technique. An illustrative example of such a system was presented. Related 
superpotentials fl29|) - fl32l) are matrices of arbitrary dimension n that can be further 
simplified with unitary transformation. 

The Hamiltonian factorization and the shape-invariance condition itself were dis- 
cussed. It was shown that an alternative representation ([5]) of operators ak and a|. 
is equivalent to a standard one (|9]). An unknown function Fk that appears in form- 
invariance condition significantly relies on the superpotential's form. If it equals to 
identical function F^ = k, the corresponding superpotential is direct sum of shifted 
one- dimensional oscillators. In other cases, it is reasonable to consider F^ to be unit 
shift, as it was shown in the discussion section. 

In spite of this, the class of obtained potentials is strongly restricted by the su- 
perpotential's form. 

In future works, a more general form of superpotential will be discussed. 
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1 Introduction 

Supersymmetric quantum mechanic presents a powerful and elegant tool for obtain- 
ing explicit solutions of quantum mechanics problems described by Schrodinger equa- 
tions [Ij. Invented by Gendenstein [2], property of discrete reparametrization of 
potentials, known as shape-invariance, helps to determine whether eigenvalues of 
Hamiltonians can be calculated by algebraic methods. 

Though shape-invariant potentials do not exhaust the full class of potentials of 
solvable Schrodinger equations \13\, it was interesting to find new integrable models. 
Previously several attempts of describing the class of shape-invariant potentials were 
made. In paper of Cooper et al [3] a wide class of scalar shape-invariant potentials 
was described. 

An attractive example of a matrix problem which admits a shape invariant su- 
persymmetric formulation was discovered by Pron'ko and Stroganov [1], who studied 
a motion of a neutral non-relativistic fermion which interacts anomalously with the 
magnetic field generated by a thin current carrying wire. The supersymmetric ap- 
proach to the Pron'ko-Stroganov problem was applied in papers ^, ^, |7j. 

Particular cases of matrix potentials were discussed in papers [8], [9], [lO], [H]. 
Matrix superpotentials appear also in some supersymmetric systems related to the 
crystalline structures in Gross-Neveu model [T7|, ^E\. Two-dimensional matrix super- 
potentials, including shape invariant, were studied in papers [I5], [I6]. In paper [11] 
Fukui considers a certain class of shape invariant potentials, which however was ad 
hoc restricted to 2 x 2 matrices with superpotentials linearly dependent on the variable 
parameter. 

Thus, in contrast to the class of scalar potentials, the class of known matrix 
potentials is presented by important but rather particular examples. The remaining 
part of the class is still undiscovered and requires further research. It seems to be 
interesting to widen the class of shape-invariant matrix potentials, because in such 
a way we will be able to describe new exactly-integrable systems of Schrodinger 
equations. 

A systematic study of the problem was carried out in our recent paper fT2] where 
we presented a complete description of irreducible matrix potentials, which include a 
term linear in variable parameter, and is proportional to the unit matrix. 

In the present paper the classification is continued and Fukui's results are gener- 
alized. All shape-invariant matrix potentials of an arbitrary dimension with super- 
potentials linearly dependent on the variable parameter were found. 

2 Shape-invariant potentials 

Let's start with a spectral problem 

Hki^ = E.^P, (1) 



where H^ is a Hamiltonian with a matrix potential, E^ and ij} are its eigenval- 
ues and eigenf unctions correspondingly, moreover -0 is n-component spinor. In the 
Schrodinger equation Hamiltonian has the form 

Hk = -^,+Vu{xl (2) 

where Vk{x) is a ra-dimensional matrix potential depending on variable x and param- 
eter k. 

Suppose that Hamiltonian accepts factorization 

Hk = alttk, (3) 

then its superpartner is defined as follows 

H^ = a^al (4) 

The most common representation of operators a^ and a|. has the form 

ak = Mx)q^ + Bk{x), a\ = -^4(^) + Bl{x), (5) 

where Ak{x)^ Bk{x) are matrices depending on x and ^^(x), B\.{x) are hermitian con- 
jugate to them. 

Substituting this representation into equation ([3]) we obtain the equation 

Hk = -4Akg^ + {BiAk - AiBk - (4^fe)')^ + BlBk - (AiBk)', (6) 

which is supposed to be Schrodinger equation of the form ([2]). It leads to the following 
conditions 

AlAk = I, 

BlAk - AiBk - (AlAkY = 0, (7) 

BlBk - {AiBk)' = Vk. 

In terms of new variable Wk{x) = Al.{x)Bk{x) this condition take the form 

Wl = Wk, 

Vk = wl -W'k. 

The same result can be obtained analogously with simpler representation of op- 
erators ak and a{. 

ak = ^ + Wk{x), ai = -— + Wk{x). (9) 



Wh{x) is called a matrix superpotential. 

As Wk{x) is hermitian, then the corresponding potential and its superpartner 
y+(x), i.e. 

v.^-?^.wl K^'-^.^n (10) 

are hermitian too. 

The goal is to find such superpotentials which generate shape-invariant Hamilto- 
nians 

Ht = Hf, + a, (11) 

where F^, Ck are scalar functions of k. In terms of superpotential, the last condition 
has the form 

Wl + Wl = W\ - W'p^ + C (12) 

It is sufficient to search for irreducible superpotentials, which means matrix Wk 
can't be transformed to block-diagonal form with a constant unitary transformation. 
Note that ii Fk = k then equation (IT2|) transforms into 

W'k = \Ck (13) 

and it follows that 

Wk = \CkX + fifc, (14) 

where Vtk is hermitian matrix depending on k which can be diagonalized. So in this 
case superpotential is completely reducible to the one-dimensional shifted oscillators. 
In case Fk ^ k it makes sense to restrict ourselves to unit shifts 

Fk = k + l, (15) 

see discussion section for details. 

In the following sections the classification of superpotentials linearly dependent 
on parameter k is presented. 



3 The determining equations 

Let us consider the superpotential of the form 

Wk = kQ + P, (16) 

where P and Q are n x n hermitian matrices dependent on x. 



Suppose that Q is not proportional to the unit matrix because in this case the 
superpotential is reducible. 

Substituting flT6|l into flT2|l and taking into account equation flT5|l we obtain the 
equation 

(2k + 1)(Q' - Q2) _ {g^ p} + 2P' = Cfc, (17) 

where Q' = '^j F" = ^ and {Q, P} = QP + PQ is an anticommutator of matrices 
Q and P. After variable separation the last equation transforms into a system of 
determining equations 

Q' = Q^ + V, (18) 

P' = \{Q.P}-li. (19) 

Cfc = (2A; + l)z/-2//, (20) 

where ^^v ^^ are arbitrary constants. For scalar values unit matrix / is omitted 
throughout the paper and written, for instance, as /i instead of [il. 

4 Solving the determining equations 

In this section a solution for determining system of arbitrary dimension n is presented. 
Let us start with equation (TT51) . 
Introducing new variable 

M = Q-Lp, (21) 

where scalar function ip is defined in the following way 

{Atan(Ax + 7), i^ = A^ > 

-Atanh(Ax + 7), i^ = -A^ < (22) 

and 7 G M is chosen so that M is not a singular matrix, we can transform equation 
(Hg) into 

M-^M'M-^ = 1 + 2ipM'\ (23) 

Note that 

(M-^y = -M-^M'M-\ 
So if we set A^ = M~^ then we obtain the following linear equation: 

N' = -I- 2ipN. (24) 



Its general solution has the form 
N = ~p{x)I + e{x)C, 



(25) 



(26) 



where p{x) and 6{x) are scalar real valued functions and C is an arbitrary constant 
matrix, 

^sin(2(Ax + 7)), z/ = A^ > 
^sinh(2(Ax + 7)), iy = -X^<0 
X + 7, z/ = 

cos^(Ax + 7), z/ = A^ > 
cosh^(A2; + 7), z/ = -A^ < 
(x + 7)2, z/ = 

Q is hermitian, so A^ is hermitian too. Moreover, if unitary matrix U is such 
that UNW is a diagonal matrix, then UQW is also diagonal. As p and 6 are scalar 
functions, then diagonalization transformation U should not depend on x. So Q can 
be diagonalized 



Q = diag{gi,...,g„} 
and equation ( [T8|) transforms into system 



(li = <lz+'^, i = l--n, 



(27) 



(2J 



which has following solutions 

gj = Atan(Aa; + 7j), i = l...n, u = X^ > 0; 
-Atanh(Ax + 7j), i = l...m 
-Acoth(Aa; + 7j), i = m + l...l , i/ 
±A, i = I + l...n 

i = l...m 
m + l...n 



Qi 



Qi 



-A^ < 0; 



(29) 



1 



0, i 



0, 



where 7j G M, i = l...n are integration constants. In cases u < and ly = matrix 
Q consists of blocks of size m,l — m+l,n — l + l and m, n — m + 1 correspondingly. 
Some of them may or may not appear, i.e. have zero size. 

So we define matrix Q up to constant unitary transformation U. Now let us 
consider linear equation flT9|) for P = {pij} which can be solved element-wise: 



If zy = A2 > 



Pa = y tan(Ax + 7i) + (pu sec{Xx + 'ji) , i = l...n 



Pij = ipijJsec{Xx + 7j) sec(Aa; + 7^), i = l...nj = l...n 

If Zy = -A2 < 



(30) 



Pa 



■J tanh(Ax + 7^) + ipu sech(Ax + 7j), i = l...m 
'J coth(Ax + 7j) + ipii csch(A2; + 7j), i = m + 1...1 



Pij 



_ ±f + ipii exp{±Xx) , 
ifij y^sech{Xx + 7j) sech(Ax + 7^) 
^Pij^y sech.{Xx + %) csch(A2; + 7^) 
ipij A/sech(Aa; + %) exp(±Ax), 
ipij ^ycsch{Xx + 7j) csch(Ax + 7^) 
ipij A/csch(Ax + 7j) exp(±Ax), 
(/?ijexp(±Ax), (*) 

'/'y> (*) 



z = / + l...n 
i = l...m,j = l...m 
i = l...m,j = 171 + 1...1 
i = l...m, J = / + l...n 
i = m + l.-.l,j = m + 1...1 
i = m+ l---l,j = I + l...n 
i = I + l...n,j = I + l...n 
i = I + l...n,j = I + l...n 



(31) 



Use formula (*) if corresponding diagonal entries qt and Qj have the same sign and 
formula (*) otherwise. 



If z/ = 



Pa 



Pij 



-/iX + Lpii, 
^(x+7i)(a:+7j)' 

Vij 

Vij^ 



i = l...m 

i = m + l...n 
i = l...m,j = l...m 
i = l...m,j = m + l...n 
i = m + l...n,j = m + l...n 



(32) 



where (fji = Tp~^ G C are integration constants. Numbers m and / in the above intervals 
correspond to the numbers defined in fl2^ . Matrix P can be further simplified with 
unitary transformation, that reduces the number of non-zero entries in similar blocks. 
Obtained superpotentials are irreducible if there are enough non-zero entries Lpji. 



5 Solving spectral problem 

Consider the Schrodinger equation 



(33) 
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where c^ vanishes with constant multiphed by unit matrix in the Hamiltonian H]^. It 
follows from shape-invariant condition f lT2|l that 



Ck = Ck+i- Ck. (34) 



Since all considered Hamiltonians are shape-invariant, equation fl33|) can be solved 
using the standard SSQM technique. An algorithm for constructing exact solutions 
of supersymmetric shape-invariant Schrodinger equations can be found in [3]. 

• Ground state '0fc(x) is proportional to the square-integrable solution of the first 
order equation 

ak^l^lix) = (^ + Wj\ ijl{x) = 0. (35) 

Function ip^ solves equation (l33l) with eigenvalues 

E'k = -Ck. (36) 

• Solution which corresponds to the n*^ excited state ipk{x) can be represented 
as 

V^^(x) = 44+1 ■ ■ ■ 4+n-iV'fc+„(a^) (37) 

and must be a square-integrable function of x. The corresponding eigenvalue is 

Ek = —Ck+n- (38) 

It follows from equations (ED]), (El]), ([2SD, dMD that 

E'^ = El + 2nfi - (n^ + 2kn)u. (39) 

At the end of the section, an example of exactly integrable problem is presented. Let 
Q and F, defined by formulas (129]). (15^. be 



_l_ 


n 


X 











IX 1 


V 


2 2x 


v^ 


\fx 


/iX 



Q=\ .^ nh ^=r ." ,r > (40) 



where ^, </? are real constants, ;U > 0. Note that, in contrast with f l32]) . the constant [i 
is changed to — /i to simplify notation. Then corresponding Schrodinger equation 
looks like 



dx^ \ir\ ^ - ^^^ ^'x^ + ^ ]\i \i 



J.3 2 



and Cfc is equal to /i. 

The square integrable ground state solution ^° = I ao ) ^^s the form 

00 = C,ci>, + C202, ...^ 

e° = Ci^l + C26 ^ ^ 

for non integer A; > 0, where 

6 = ^x^^'e-'^H'^iK 0,-2 -k,'^, f x)+ (43) 



and 



2^^ ^ -'-'OV'",", -- •", ^, 2 



xh-^HB{-k,0,-2-k,^,^x), 



^' 2 
6 = f xe-^if^(-fc, 0, -2 - fc, ^, f x)+ (44) 

^e-^//M-A:,0,-2-fc,^,fx). 

Here iif^ denotes the Heun Biconfluent function. For integer k solutions ( H3|l and 
(jH]) are identical and second linearly independent solution of equation fl35l) should 
be used instead of fH4|) . It is however cumbersome and is omitted in the paper. 
Note that constants Ci and C2 should be chosen so, that 

o 

(i^y^j^dx = 1. (45) 

The spectrum of the problem is described by formula 

E" = -(2n + l)/i. (46) 

6 Discussion 

In this section the possibility of transforming unknown function F^ in the shape- 
invariance condition flTTl) to unit shift is discussed. 

During the computation it was assumed that Fk = k + 1, let us show it is a 
reasonable restriction. Under invertible transformation of variables 

k -^ a{k) (47) 

the function Fj. = F{k) changes by a similar transformation 

F{k) -^ aFa-\k) = a{F{a-\k))). (48) 



Searching for such transformation, that would change function F{k) to unit shift, we 
get the equation 

F{a-\k)) = a~\k + l). (49) 

The above equation is known as Abel functional equation. The results concerning 
the solution of this equation were obtained in papers [H], [20], [21], [22]. 

Let X be M or M+. It is proved that 

(CI) if F : X ^ W is an injective function such that for every compact set K <Z X 
there exists p G N such that Vn, m G N°, \n — m\ > p : 

F'\K) n F'^iK) = 

then there exists a solution to the Abel functional equation. 

So if the above condition is fulfilled, F^ can be transformed into unit shift. 

7 Conclusion 

Restricting the class of superpotentials to a simple form f lTB]) . it was possible to find 
new exactly solvable systems of Schrodinger equations. As corresponding Hamiltoni- 
ans satisfy shape-invariant condition (IT2|) . the systems can be solved using a standard 
SSQM technique. An illustrative example of such a system was presented. Related 
superpotentials (!29|) - (!32l) are matrices of arbitrary dimension n that can be further 
simplified with unitary transformation. 

The Hamiltonian factorization and the shape-invariance condition itself were dis- 
cussed. It was shown that an alternative representation ([5]) of operators ak and a|. 
is equivalent to a standard one (|9]). An unknown function Fk that appears in form- 
invariance condition significantly relies on the superpotential's form. If it equals to 
identical function F^ = k, the corresponding superpotential is direct sum of shifted 
one-dimensional oscillators. In other cases, it is reasonable to consider F^ to be unit 
shift, as it was shown in the discussion section. 

In spite of this, the class of obtained potentials is strongly restricted by the su- 
perpotential's form. 

In future works, a more general form of superpotential will be discussed. 
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